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ABSTRACT 

We provide a general analysis on the properties of emitting material of some rapidly evolving and 
luminous transients discovered recently with the Pan-STARRSl Medium Deep Survey. It is found 
that these transients are probably produced by a low-mass non-relativistic outflow that is continuously 
powered by a newly born, rapidly spinning, and highly magnetized neutron star. Such a system could 
originate from an accretion-induced collapse of a white dwarf or a merger of a neutron star-neutron 
star binary. Therefore, observations to these transients would be helpful for constraining white dwarf 
and neutron star physics and/or for searching and identifying gravitational wave signals from the 
mergers. 

Subject headings: gamma-ray burst: general — stars: neutron — supernovae: general 


1. INTRODUCTION 

Neutron stars (NSs) are usually born from core col¬ 
lapse of a massive star, during which a luminous super¬ 
nova can be driven by the intense neutrino emission from 
the proto-NS by expelling and heating the stellar enve¬ 
lope. Besides this most popular origination, there could 
be other two ways leading to a NS formation. On one 
hand, it is believed that, if electron captures in the core 
can take place more quickly than the nuclear burning, 
an accreting white dwarf (WD) could collapse into a NS 
by loosing pressure support rather than lead to a disrup¬ 
tion event (e.g. Type la supernova; Canal & Schatzman 
1976; Nomoto & Kondo 1991). On the other hand, re¬ 
cently it is often suggested that a massive NS could form 
and survive from a merger of a NS-NS binary (Dai et 
al. 2006; Fan & Xu 2006; Metzger et al. 2008; Buc- 
ciantini et al. 2012; Rowlinson et al. 2010, 2013), during 
which a highly-beamed relativistic jet can be launched to 
produce a short-duration gamma-ray burst (GRB; Nakar 
2007; Berger 2014). The clues to the post-merger NSs 
have emerged more and more from the observations of 
extended emission, X-ray flares, and plateau afterglows 
of some short GRBs. 

During accretion-induced collapses of WDs and NS- 
NS mergers, a low-mass non-relativistic outflow can 
be ejected nearly isotropically due to the dynamical 
centrifugal force and the wind blown from a remnant 
NS+disk system. In such an outflow, some radioactive 
heavy elements are expected to be synthesized. Then 
the decay of the elements would generate a faint optical- 
infrared transient with a duration of several days and 
a luminosity of < 10 41 erg s _1 (Darbha et al. 2010; Li 
& Paczynski 1998; Kulkarni 2005; Rosswog 2005; Met¬ 
zger et al. 2010; Roberts et al. 2011). In the NS-NS 
merger case, such transients are now known as kilono- 
vae or macronovae, which have become one of the recent 
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highlights of time domain astronomy since from the ob¬ 
servation of the infrared bump after GRB 130603B (Tan- 
vir et al. 2013; Berger et al. 2013). 

However, as implied by the X-ray and optical af¬ 
terglows of GRB 130603B, the GRB jet and thus the 
isotropic outflow could both be primarily powered by the 
spin-down of a newly born NS rather than radioactivi¬ 
ties (Fan et al. 2013), although its infrared luminosity is 
still too ambiguous to achieve a definite conclusion. In 
any case, this observation encouraged us optimistically to 
consider that the luminosity of the transients during WD 
collapses and NS-NS mergers could be enhanced signif¬ 
icantly by a spinning-down post-collapse/merger NS to 
be comparable to or even exceed those of ordinary su¬ 
pernova, but still with an obviously short duration due 
to the low mass of the outflow (Yu et al. 2013; Met¬ 
zger & Piro 2014). Reasonably, a new class of rapidly 
evolving and luminous transients were recently reported 
by Drout et al. (2014) with the Pan-STARRSl Medium 
Deep Survey (PS1-MDS), which could therefore be the 
first smoking-gun observational signatures for the NS- 
powered transient events during WD collapses and/or 
NS-NS mergers. 

2. THE PS1-MDS TRANSIENTS 

The PS1-MDS transients were specifically identified by 
following criteria: the transient must rise by >1.5 mag in 
the 9 days immediately before a peak and decline by >1.5 
mag in ^25 days after the peak. Such rapid evolution of 
the transients makes them easily evade from previous su¬ 
pernova surveys until the PS1-MDS which is featured in 
its rapid cadence and multiple band coverage to a signifi¬ 
cant depth (~24 mag). For 10 of the new-discovered PS1- 
MDS transients, the spectra of their underlying hosts 
were obtained so that the cosmological redshifts of the 
transients can be determined to range from z = 0.074 
(PSl-lOah) to 2 = 0.646 (PSl-llbbq). Meanwhile, the 
evolution of the spectral energy distributions of the tran¬ 
sients selves was also achieved by interpolating their light 
curves in five broadband filters (gpiTpiipi^pi i/pi). Then, 
the effective black-body temperatures of the transients 
(except for PSl-12bb) were revealed to evolve around 
20,000 K, which indicates an ultraviolet emission, near 
the peak emission time to 7000 K at later times. The 


2 


pseudo-bolometric luminosities of the transients were 
correspondingly calculated to have a peak value with the 
order of L p ~ 10 43 erg s” 1 . Consequently, the radius of 
emission region at the peak time can be constrained to be 

R p = (L p /47rcrTp) 1/2 - 3.0 x 10 14 cm L^ 3 T~l 3 , where 
a is the Stefan-Boltzmann constant and, hereafter, the 
conventional notation Q x = Q/10 x is adopted in the cgs 
units. 

The above radius, luminosity, and temperature of the 
PS1-MDS transients indicate that they are probably as¬ 
sociated with stellar-size energetic “explosions”. Ac¬ 
cording to the selection criteria of Drout et al. (2014), 
the peak times t p of the transient light curves can be 
taken to be from a few to ~ 10 days after the explo¬ 
sions. Such short peak times make the transients dis¬ 
tinct from the overwhelming majority of ordinary su¬ 
pernovae, which usually reach their peaks after several 
weeks or a few months from the explosions. More pre¬ 
cisely and specifically, the PS1-MDS transients with an 
intermediate luminosity < 10 43 erg s _1 (i.e. PSl-lOah, 
PSl-lObjp, PSl-12bb, PSl-12brf, PSl- 13dwm, and PS1- 
13ess) still exhibit some similarities with the previously- 
reported rare supernova events such as SN 2002bj (Poz- 
nanski et al. 2010), SN 2010X (Kasliwal et al. 2010), 
and SN 2005ek (Drout et al. 2013). However, the com¬ 
pletely new luminosity-timescale region opened by the 
high-luminosity PS1-MDS transients (i.e. PS1- llqr, 
PSl-llbbq, PSl-12bv, and PSl-13duy), which has never 
been probed before, would lead us to consider a com¬ 
pletely new branch of stellar explosions. 

Regardless of the specific nature of the progenitors of 
the PS1-MDS transients, a general order-of-magnitude 
analysis can be made with respect to a cluster of mate¬ 
rial of mass M and internal energy £i nt expanding due to 
its own internal pressure. On one hand, the peak of the 
material emission generally appears at the time within 
which photons inside the material can transfer from the 
core to the outmost surface. By considering of the ran¬ 
dom walk of the photons, the peak emission time can be 
expressed by t v ~ (R p / A) 2 (A/c) = 3kM/4ttR p c , where 
A = («p) -1 is the average free path, n the opacity, and 
p the density. Then the mass and the kinetic energy of 
the material can be estimated by 


= °' 04 M ° *= 0 . 3 ^ 43 ^ 4 . 3 ^ 6 , ( 1 ) 


and 


E Kp = -Mvl ~ 3.3 x 10 48 erg (2) 


respectively, where the expanding velocity is given by 


v p 


^ = 3.0 x 10 8 cm s -1 


T 1/2 rji — 2 J.— 1 
-^p,43- t p,4.3 6 p,6- 


(3) 


On the other hand, the emission luminosity can be ap- 
proximatively connected with the internal energy by 
L ~ AnREintc/ (3ftM) by following the diffusion equa¬ 
tion / = (c/3Kp)(dui nt /dr), where / = L/(47ri? 2 ) is the 
emission flux, Ui„t = -Ednt/V the internal energy density, 
and V = 17ri? 3 the volume. Then we get 

_ 3 kML v aq r . . 

A’int.p ~ ~^R~c = = X er § (^) 


TABLE 1 

Parameters for fittings to the light curves 


Events 

M/ 10- 2 M q 

£Asd,i/10 44 erg s" 1 

^■sd / day 

PSl - llbbq 

2.0 

4.5 

2.0 

PS1 - llqr 

1.5 

1.5 

5.0 

PSl - 12bv 

1.5 

0.65 

14 


Together with the equation 4irR 2 p/M = dv/dt ~ v/t 
describing the adiabatic acceleration of the material, the 
expanding velocity at the peak time can also be calcu¬ 
lated by 


4nR 2 p p t p „ , 

-^-= 4.5 x 10 8 cm s 1 k 

M 


-0.3^p,4.3*p,6> ( 5 ) 


where the pressure p = Aui nt is adopted. By combining 
Equations © and we can obtain a rough constraint 
on the opacity asu~ 0.7 cm 2 g _1 L^J 3 T~^ 3 t~^ 7 which 
indicates an opacity primarily contributed by electron 
scattering. The above constraint also implies that the 
basic properties (e.g. mass, expanding velocity, and in¬ 
ternal energy) of adiabatically expanding material can 
be certainly determined with the peak values of its bolo- 
metric light curve. 


3. ENERGY ENGINE 

The mass of the emitting material, which is dramat¬ 
ically smaller than those of ordinary supernova ejecta, 
suggests that the PS1-MDS transients probably involve 
an explosion of a system governed by an extremely- 
stripped star or, more probably, stellar-size compact ob¬ 
jects. The material ejected during the explosion could 
be heated by the energy engine “impulsively” or by some 
other energy sources after the explosion. 

First of all, as a conventional consideration, the emit¬ 
ting material could be impulsively heated by a shock 
and therefore the emission luminosity is determined by 
the shock breakout radius. However, owing to a limit 
on the radius of the largest red supergiant (Levesque et 
al. 2009), only the PS1-MDS transients with an inter¬ 
mediate luminosity < 10 43 erg s _1 can possibly be at¬ 
tributed to a shock breakout from a normal star. For the 
high-luminosity ones (i.e. PS1- llqr, PSl-llbbq, PS1- 
12bv, and PSl-13duy), Drout et al. (2014) suggested 
that these transients could be associated with a shock 
breakout within a dense stellar wind. In this case, how¬ 
ever, an unexpected high mass-loss rate and an ad hoc 
mass-loss history of the progenitor need to be invoked. 
Then, Kashiyama & Quataert (2015) proposed that the 
impulsive heating could be accomplished in a very hot 
super-Eddington fall-back accretion disk surrounding a 
newly born black hole, where a low-mass outflow can be 
driven by the disk radiation pressure. Such a model can 
work if the collapsing progenitor star is elaborately de¬ 
signed to eject a just-right mass of material so that a 
successful supernova can be avoided but an appropriate 
fail-back accretion disk can still be maintained. 

In contrast, a post-explosion energy supply could 
be a more natural and attractive choice, which can 
be provided by radioactive decay of heavy elements 
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Fig. 1.— Comparison between the model-predicted (solid lines) and the observational pseudo-bolometric luminosity, radius, and tem¬ 
perature of the high-luminosity PS1-MDS transients (solid circles). The source PSl-13duy is excluded because only two data points were 
provided. 


and/or a long-lasting energy engine. In the radioactiv¬ 
ity scenario, the mass of heavy elements represented by 
56 Ni is required to be Mm ~ mNi-Eint, P /(?Ni + <?Co) ~ 
0.05 M 0 L Pi 43 t Pi 6 , where WNi is the mass of a 56 Ni atom 
and (?Ni =1-78 MeV and qc 0 =3.505 MeV are the energy 
release for the decay of per 56 Ni and per subsequent 56 Co, 
respectively. As shown, for a sufficiently high luminosity, 
the required nickel mass could unrealistically exceed the 
total mass of the emitting material. Even for a relatively 
low luminosity, an unacceptable high nickel fraction is 
also predicted, because the nickel mass is only slightly 
lower than the total mass. Therefore, the radioactivity 
must not be the primary energy source for the PS1-MDS 
transients and alternatively a long-lasting energy engine 
becomes necessary and crucial. The long activity of the 
engine further requires that the post-explosion compact 
object should be a NS rather than a black hole, as previ¬ 
ously suggested for GRBs (Dai & Lu 1998a,b; Zhang & 
Meszaros 2001) and superluminous supernovae (Kasen 
& Bildsten 2010). The energy that a NS can provide 
to power the explosion outflow mainly comes from its 
spin-down with a luminosity L s d(t) = L s d,i (1 + t/t s d) 2 • 
Then by taking appropriate values of the parameters 
M, L s d,i and t s( j as listed in Table 1, we employ a 
semi-analytical light curve model (see Appendix for de¬ 
tails) to account for the pseudo-bolometric light curves of 
the high-luminosity PS1-MDS transients (i.e., PSl-llqr, 
PSl-llbbq, and PSl-12bv) as well as the evolution of 
radius and temperature of the transients, as presented 
in Figure [D In our fittings, a typical constant opac¬ 
ity k = 0.5cm 2 g _1 is artificially adopted to reduce the 
parameter freedom and an uncertain parameter £ is in¬ 
troduced to represent the energy injection efficiency. 

Being represented by the braking mechanism due to 
magnetic dipole radiation, the spin-down luminosity 
and timescale of a NS can be expressed by A s d,i = 


TABLE 2 

Neutron star parameters 


Events 

Pi/ms 

S/10 14 G 

PSl - llbbq 

me /2 

17^/ 2 

PS1 - llqr 

18£P 2 

12£ 1//2 

PSl - 12bv 

165 1 / 2 

6.5? 1 / 2 


10 7 erg s -1 S 2 Pj -4 and t sd = 2 x 10 39 s P _2 P 2 , where 
B and P; are the magnetic field strength and initial spin 
period, respectively. Then, from the fitting parameters, 
the values of B and P, can be derived as functions of 
the parameter £, as presented in Table 2. It is shown 
that the involved NSs should be highly magnetized and 
spin rapidly with a millisecond period, in spite of the un¬ 
certainty of the energy injection efficiency. In short, the 
PS1-MDS transients (at least the high-luminosity ones) 
are probably produced by a low-mass outflow that is con¬ 
tinuously powered by a millisecond magnetar. Specif¬ 
ically, such a newly born NS could originate from an 
accretion-induced collapse of a WD or a merger of a NS- 
NS binary. During these processes, a strong magnetic 
field can be naturally generated and amplified via flux 
freezing and dynamo action (Duncan & Thompson 1992; 
Cheng & Yu 2014). On the contrary, for an extremely- 
stripped star, it is nearly impossible to produce a mil¬ 
lisecond magnetar through its collapse. 

4. DISCUSSIONS 

In view of the high similarity in the post-explosion con¬ 
figurations between the WD collapses and NS-NS merg¬ 
ers, it is not easy to clearly discriminate these two models 
by the present observations of the PS1-MDS transients. 

Firstly, the derived spin periods of ~ 16 — 18 ms for 
£ ~ 1 seems only prohibited by WD collapses but con- 
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flict with NS-NS mergers, because the high initial angu¬ 
lar momentum of the binary would give a limiting spin 
period of P; ~ 1 ms to the post-merger NSs. However, 
in fact, it is widely believed that such a near-Keplerian 
spinning NS would lose its rotational energy immediately 
by radiating gravitational waves (Dall’Osso et al. 2015), 
amplifying magnetic fields (Cheng & Yu 2014), and pro¬ 
ducing short GRBs and subsequent extended emission 
(Metzger et al. 2008; Bucciantini et al. 2012). As a 
result, the initial period for magnetic dipole radiation 
would become much longer than ~ 1 ms, e.g., P\ = 5 
ms (e.g. Cheng & Yu 2014). For such a reference value, 
we can constrain the energy injection efficiency to be 
£ ~ 10%. This might be caused by that (1) the major¬ 
ity of the spin-down luminosity is concentrated into the 
GRB direction to contribute a plateau afterglow emission 
(Yu et al. 2010), which however probably deviates from 
the light of sight and, perhaps, (2) the spin-down lumi¬ 
nosity on the light of sight could also be partly released 
through X-ray emission by the NS wind itself (Dai 2004; 
Metzger & Piro 2014). Therefore, such wind X-ray emis¬ 
sion associated with the UV-optical transients, if exists, 
could provide a testable signal for our model. 

Secondly, the NS-NS merger model is seemingly chal¬ 
lenged by other two issues. (1) As pointed out by Drout 
et al. (2014), the explosion site offsets of all tran¬ 
sients with respect to their host galaxies are statistically 
smaller than those expected for NS-NS mergers. (2) As 
calculated by Kasen et al. (2013), the opacity of the 
neutron-rich outflow could be on the order of magnitude 
of 10 — 100cm 2 g -1 due to the synthesization of lan¬ 
thanides, which is much higher than the opacity inferred 
here. However, on one hand, if we individually pick up 
the offsets of the high-luminosity transients, we will find 
that most of them are actually on the high side of the 
offset distribution. On the other hand, the opacity of the 
outflow actually can be reduced by two effects: (1) the 
lanthanide synthesization in the disk wind component 
of the outflow is blocked due to the neutrino irradiation 
from the post-merger NS (Metzger & Fernandez 2014) 
and (2) the lanthanides in the dynamical component of 
the outflow could be ionized by the X-ray emission from 
the NS wind (Metzger & Piro 2014). 

Finally, a potential discrimination between the two 
models may arise from the explanation of the event rate 
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of PS1-MDS transients, i.e., roughly, ~ 3200Gpc _3 yr _1 
corresponding to the four high-luminosity transients. As 
a comparison, the rate of Type la supernovae was found 
to (3.01 ± 0.062) x 10 4 Gpc~ 3 yr _1 from the Lick Obser¬ 
vatory Supernova Search (Li et al. 2011) and the rate of 
NS-NS mergers is also widely considered to be on the or¬ 
der of ~ 1000Gpc _3 yr _1 (Kalogera et al. 2004). There¬ 
fore, optimistically, the event rate of PS1-MDS transients 
could in principle be accounted for in both the WD col¬ 
lapse and NS-NS merger models, as long as a remark¬ 
able fraction of these processes can lead to a NS forma¬ 
tion. The rationality of such fractions could be tested by 
some further model implications, e.g., the cosmic frac¬ 
tion of neutron-rich elements, the detectability of as¬ 
sociated gravitational wave signals by Advanced LIGO 
and Virgo (Abadie et al. 2010; Nissanke et al. 2013), 
and the detectability of subsequent radio transients aris¬ 
ing from the interaction between the outflow and circum 
medium (Metzger et al. 2015), etc. For this purpose, 
some detailed calculations and simulations need to be 
implemented with the constraints derived in this Letter. 

5. CONCLUSION 

The luminosity-timescale region opened by the PS1- 
MDS transients, which although could be partly over¬ 
lapped by the shock breakout emission of some super¬ 
novae, leads us to get a glimpse of a completely new 
branch of stellar-size explosion phenomena. The mass 
and energy requirements on the emitting material for 
the transients robustly indicate a post-explosion system 
consisting of a newly born millisecond magnetar and a 
low-mass non-relativistic isotropic outflow. Such a sys¬ 
tem could originate from accretion-induced collapses of 
WDs or mergers of NS-NS binaries. Therefore, future 
observations to these transients would play a vital role 
in constraining WD and NS physics and/or in searching 
and identifying gravitational wave radiation. 
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APPENDIX 

THE LIGHT CURVE MODEL 


Following Arnett (1980) and Kasen & Bildsten (2010), a semi-analytical light curve model for NS-powered merg- 
ernovae was provided by Yu et al. (2013), where relativistic effects were taken into account. Here, in view of the 
non-relativistic velocity inferred from the PS1-MDS transients, we alternatively adopt a simplified non-relativistic ver¬ 
sion of the model, which was previously described in detail by Kasen & Bildsten (2010) for superluminous supernovae. 
In this case, the dynamical evolution of an outflow can be determined by 


dv 4:irR 2 p 

dt M 


and 



(A2) 


where the pressure p = i(ui nt — u e ) with zij nt being the internal energy density and u e being the energy density 
releasing by emission. With an energy supply rate £L S d, the evolution of the outflow internal energy is determined by 
the energy conservation law as 


dE lnt 

dt 


£L sd - L -p 


dV 
dt ’ 


(A3) 


where L is the emission luminosity and the work —pdV = —47T R 2 pv represents the energy loss due to the expansion 
of the outflow. The emission luminosity can further be related to the internal energy by L = E[ nt c/(rR) for t < t T 
and L = E- lnt c/R for t > t Tl where r = 3kM/AttR 2 is the optical depth and the time t T represents the time for the 
optical thick-thin transition. Then, by setting the parameters k, M, £L sdi j, and t sd , a light curve can be solved from 
the above equations. Three typical types of light curves are presented in Figure [21 
For an analytical understanding of the light curve profiles, we can approximate the spin-down power of the NS by 
a broken-power-law function as 


L s d —£A d; i, for t <C fsd> (A4) 

= ££d,i(iAsd) -2 , for t > t sd . (A5) 

Meanwhile, by introducing the photon diffusion time of td = (2>kM/A nvc) 1 ^ 2 with a constant velocity, Equation (IA3I) 
can be rewritten to (Kasen & Bildsten 2010) 


dE\ n R 

dt 


E. 


'supi^d 


l tl - E\ 


intC 


(A6) 


for t <t T . The shape of light curves mainly depends on the relation between the timescales t d and t sd . Let us consider 
the case of £ sd < f d as an example. Firstly, for t < f sd , the energy supply is constant and the solution of Eq. (IA6I) can 
be derived to 


L = 





(A7) 


Secondly, for f sd < t < £ d , by omitting the minor important radiation effect and simplifying Eq. (IA6I) to d(Ei nt t)/dt r; 
CL d ,it 2 d /t, we can get 


L 



1 _ e -*sd/2td 


+ 


e^dA; 


sd 



^-sd 


£-kd,i Ad 


t 


2 

d 



(A8) 
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Fig. 2.— Three typical types of light curves, where the approximative temporal behaviors are labeled for different evolution phases. 
Thirdly, for t > t d, the radiation term becomes dominative in Eq. (IA6I) . i.e., d(E- ln tt) / dt « — -Ei nt t 2 /t 2 , which gives 


L 


£LdA (l - e 


e^dhr^L 

^sd. 




2 

sd 






(A9) 


Finally, as the emission luminosity deceases to approach the energy supply luminosity, the succeeding evolution of 
the emission luminosity should track the evolution of the energy supply (i.e., L oc t ~ 2 ), because after t<j the freshly 
deposited internal energy can escape from the outflow nearly immediately. The above treatment can be easily extended 
to the other two cases of t s( j ~ td and t s( j > td . Then a completed solution to Eq. (IA6I) can be summarized as follows: 

Case I: For t s d < td 



L oc | 

r t 2 , for t < tsd, 

lnt, for t s d < t < td, 
e - * 2 , for td < t < t c , 
k t -2 , for t > t c , 

(A10) 

Case II: For t s d ~ td 

T „ f t2 , for t < t sd , 

\t~ 2 , for t> t sd , 

(All) 

Case III: For t s a > td 

L oc < 

f t 2 , for t < t d , 

t°, for t d < t < t sd , 

[ t -2 , for t > tsd, 

(A12) 


where t c is defined as the time at which L —>• £L S( j. The above three different solutions just correspond to the three 
types of light curves presented in Figure [2] In all cases, the peak luminosity of the emission appears at the diffusion 
time td- The three high-luminosity PS1-MDS transients are all fitted by the model with t s d ~ td in the text. 















